arXiv:hep-th/0512221v3 5 Mar 2006 


One loop boundary effects: techniques and 

applications. * 

Valery N.Marachevsky ^ 

V. A. Fock Institute of Physics, St. Petersburg University, 
198504 St. Petersburg, Russia 

February 1, 2008 


Abstract 

A pedagogical introduction to the heat kernel technique, zeta func¬ 
tion and Casimir effect is presented. Several applications are consid¬ 
ered. First we derive the high temperature asymptotics of the free 
energy for boson fields in terms of the heat kernel expansion and zeta 
function. Another application is chiral anomaly for local (MIT bag) 
boundary conditions. Then we rederive the Casimir energies for per¬ 
fectly conducting rectangular cavities using a new technique. The new 
results for the attractive Casimir force acting on each of the two per¬ 
fectly conducting plates inside an infinite perfectly conducting waveg¬ 
uide of the same cross section as the plates are presented at zero and 
finite temperatures. 


1 Introduction 

The main problem of the quantum field theory with boundaries is its renor¬ 
malization and physical meaning of the results obtained. Divergences that 
appear in quantum field theory make the problems on manifolds with bound¬ 
aries more complicated than in infinite space. 

In the presence of boundaries or singularities the heat kernel technique 
is an effective tool for the analysis of the one loop effects (see reviews P, 
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j2]). Different applications of the heat kernel expansion exist. The heat 
kernel techniqne seems to be the easiest way for the calcnlation of qnantum 
anomalies, calculation of effective actions based on finite-mode regularization 
and analysis of divergences in quantum field theory. 

Chiral anomaly, which was discovered more than 35 years ago j3j, still 
plays an important role in physics. On smooth manifolds without boundaries 
many successful approaches to the anomalies exist The heat kernel 

approach to the anomalies is essentially equivalent to the Fujikawa approach 
jlj and to the calculations based on the finite-mode regularization jH], but 
it can be more easily extended to complicated geometries. The local chiral 
anomaly in the case of non-trivial boundary conditions (MIT bag boundary 
conditions) has been calculated only recently p. 

Casimir effect PI is a macroscopic quantum effect. Briefly speaking, if we 
impose classical boundary conditions on a quantum field on some boundary 
surface than we get the Casimir effect. There are several different physical 
situations that should be distinguished in the Casimir effect. 

Suppose there are two spatially separated dielectrics, then in a dilute 
limit (e —>■ 1) the Casimir energy of this system is equal to the energy of 
pairwise interactions between dipoles of these two dielectrics via a Casimir- 
Polder retarded potential m- For a general case of separated dielectrics the 
Casimir energy can be calculated as in PI or pn] (for a recent discussion of 
these issues see [H] and a review PI. new possible experiments in uni ). 

A different situation takes place when there is a dilute dielectric ball or 
any other simply connected dielectric under study (see a review [TJj for a 
discussion of related subjects and methods used). As it was pointed out in 
PI and then discussed in detail in PI. microscopic interatomic distances 
should be taken into account to calculate the Casimir energy of a dilute 
dielectric ball. The average interatomic distance A serves as an effective 
physical cut off for simply connected dielectrics. 

The limit of a perfect conductivity (e —> -|-cxd) is opposite to a dilute case. 
This is the strong coupling limit of the theory. Any results obtained in this 
limit are nonperturbative ones. 

The Casimir energy for a perfectly conducting rectangular cavity was hrst 
calculated in PI using exponential regularization. Later it was derived by 
some other methods (see references and numerical analysis in PI. also a re¬ 
view PI)- In the present paper we derive the Casimir energy for rectangular 
cavities at zero temperature by a new method described in Sec.3. By use 
of this method we could rewrite the Casimir energy for rectangular cavities 
in the form that makes transparent its geometric interpretation. Also this 
method yields new exact results for the Casimir force acting on two or more 
perfectly conducting plates of an arbitrary cross section inside an inhnite 
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perfectly conducting waveguide of the same cross section. 

The paper is organized as follows. In Sec.2 we give an introduction to the 
formalism of the heat kernel and heat kernel expansion. Also we introduce a 
^-function pa and calculate the one loop effective action in terms of (- func¬ 
tion. Then we consider two examples. First we derive the high temperature 
expansion of the free energy for boson fields [21] in terms of the heat kernel 
expansion and (^-function. Then we derive a chiral anomaly in four dimen¬ 
sions for an euclidean version of the MIT bag boundary conditions jUj. Sec.3 
is devoted to the Casimir effect for perfectly conducting cavities. In Sec.3.1 
we introduce a regularization and a convenient new method of calculations 
using an example of two perfectly conducting parallel plates, then apply it 
to more complicated rectangular geometries. We present the Casimir energy 
of the cavity in the form (jnj). Then we discuss an argument principle and 
C-functional regularization for the cavity. In Sec.3.6 we rewrite the Casimir 
energy of the cavity in terms of geometric optics (Unni). In Sec.3.7 we describe 
a possible experiment and derive the formula for the attractive force acting 
on each of the two parallel plates inside an infinite rectangular waveguide 
with the same cross section. Also we present a generalization of this result 
for the case of the two parallel perfectly conducting plates of an arbitrary 
cross section inside an infinite perfectly conducting waveguide with the same 
cross section at zero and hnite temperatures. 

2 Spectral techniques 

2.1 Heat kernel 

Consider a second order elliptic partial differential operator L of Laplace type 
on an n-dimensional Riemannian manifold. Any operator of this type can be 
expanded locally as 

L = + 6 ), ( 1 ) 

where a and b are some matrix valued functions and is the inverse metric 
tensor on the manifold. For a flat space g'^‘' = 6^'^. 

The heat kernel can be defined as follows: 

K{t]X]y]L) = {x\e-x.Y>{-tL)\y) = ^ 0j^(a;)0A(2/) exp(-fA), (2) 

A 

where 0 a is an eigenfunction of the operator L with the eigenvalue A. 

It satisfies the heat equation 

{dt +L^)K{t-,x-,y,L) = Q (3) 
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with an initial condition 


K{Q-x-y-L) = 5{x,y). (4) 

If we consider the helds in a hnite volume then it is necessary to specify 
boundary conditions. Different choices are possible. In section 3.1 we will 
consider the case of periodic boundary conditions on imaginary time coor¬ 
dinate, which are specihc for boson helds. In section 3.2 we will study bag 
boundary conditions imposed on fermion helds. If the normal to the bound¬ 
ary component of the fermion current vanishes at the boundary, one 

can impose bag boundary conditions, a particular case of mixed boundary 
conditions. We assume given two complementary projectors n±, n_-|-n+ = I 
acting on a multi component held (the eigenfunction of the operator L) at 
each point of the boundary and dehne mixed boundary conditions by the 
relations 

n_t/’|aM = o, (v„-f5)= 0, (5) 

where S' is a matrix valued function on the boundary. In other words, the 
components n_'^ satisfy Dirichlet boundary conditions, and 11+'^ satisfy 
Robin (modihed Neumann) ones. 

It is convenient to dehne 


X = n+ - n_. (6) 

Let {cj}, j = 1,..., n be a local orthonormal frame for the tangent space 
to the manifold and let on the boundary be an inward pointing normal 
vector. 

The extrinsic curvature is dehned by the equation 

Lat = K, , (7) 

where T is the Christohel symbol. For example, on the unit sphere 
which bounds the unit ball in i?" the extrinsic curvature is Lab = Sab- 

Curved space ohers no complications in our approach compared to the 
hat case. Let Rfj,upa be the Riemann tensor, and let = R"^pua be the 
Ricci tensor. With our sign convention the scalar curvature R = R^ is +2 
on the unit sphere S'^. In hat space the Riemann and Ricci tensors are equal 
to zero. 

One can always introduce a connection Up and another matrix valued 
function E so that L takes the form: 

L = -{g^''VpV, + E) ( 8 ) 
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Here is a sum of covariant Riemannian derivative with respect to metric 
g^i, and connection o;^. One can, of course, express E and oj in terms of a^, 
b and 


(9) 

E = b- g^''{dyUJ^ + - ujpT^^^) (10) 

For the future use we introduce also the field strength for ui\ 

du^^ F • (H) 

The connection will be used to construct covariant derivatives. The 
subscript ; g,.. .ua will be used to denote repeated covariant derivatives 
with the connection u and the Christoffel connection on M. The subscript 
: a.. .be will denote repeated covariant derivatives containing u and the 
Christoffel connection on the boundary. Difference between these two co¬ 
variant derivatives is measured by the extrinsic curvature ©• For example, 

E-^ab E-af, LcibE-pi- 

Let US define an integrated heat kernel for a hermitian operator L by the 
equation: 

K{Q, L,t) := Tr (Q exTp{—tL)) = f (Ex^/gii {Q{x)K{t-,x-,x-,L)) , (12) 

J M 

where Q{x) is an hermitian matrix valued function, tr here is over matrix 
indices. For the boundary conditions we consider in this paper there exists 
an asymptotic expansion [2B] as f —> 0: 


R(Q, L,t)e,J2 MQ, . (13) 

k=0 

According to the general theory [2^] the coefficients ak{Q,L) are locally 
computable. This means that each ak{Q,L) can be represented as a sum of 
volume and boundary integrals of local invariants constructed from Q, D, E, 
the curvature tensor, and their derivatives. Boundary invariants may also 
include S, Lab and x- Total mass dimension of such invariants should be k 
for the volume terms and k — 1 for the boundary ones. 

At the moment several coefficients of the expansion m are known for 
the case of mixed boundary conditions © and matrix valued function Q (see 
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P for details of derivation; the formula (jSII) for 04 was derived in P with 
additional restrictions Lab = 0 and S' = 0) : 


ao(Q, L) 

= {A7r)--R [ d-x^triQ). 

J M 

(14) 

ai{Q,L) 

= h4T,yin-i)/2 f d^-^xVhtiixQ). 

(15) 


4 JdM 


(^2{Qi L) 

= ^(47r)-"/2 1 y d-x^tr (6QE + QR) 



+ f d'^~^X^/htT {2QLaa + l‘2QS + 3xQ-n)\ ■ 

(16) 


JdM ) 


L) 

= ^(47r)-(”-i)/2 f d'^-^x\/htT{Q{-24:E + 2AxEx 
384 



+48x^ + 48 Ex - 12X:aX:a + 12X:aa “ 6X:aX:aX + 

-f-SX-Ranan + 192S^ + 96LaaS + (3 + 10x)-haa-hbf, 

+ (6 — Ax)LabLab) + Q;n(96S' + 1925'^) + 2AxQ-nn}- (17) 

For a scalar function Q and mixed boundary conditions the coefficients 04 
and 05 were already derived m 

2.2 (^-function 

Zeta function of a positive operator L is dehned by 

ClW = Zv' 

A 

where the sum is over all eigenvalues of the operator L. The zeta function is 
related to the heat kernel by the transformation 

-I ^ + CXD 

a{s) = -— dtR-^K{I,L,t). (19) 

r(s) Jo 

Residues at the poles of the zeta function are related to the coefficients of 
the heat kernel expansion: 

ak{I, L) = ReSs=(n-fc)/2(r(s)CL(s)). (20) 

Here / is a unit matrix with a dimension of the matrix functions in 
(CD .From (unD it follows that 


ari(7, L) = Cl(0) . 


( 21 ) 
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In Euclidean four dimensional space the zero temperature one-loop path 
integral over the boson helds cj) = Cx(j)x can be evaluated as follows (up 
to a normalization factor): 

Z = f - n / f^dCxe-^^" ~ _ (22) 

J ^ J 

Here we introduced the constant /x with a dimension of mass in order to keep 
a proper dimension of the measure in the functional integral. Cl( 0) can be 
thought of as a number of eigenvalues of the operator L. For the operator 
L in the form (HD the number of eigenvalues is inhnite, so Cl( 0) yields a 
regularized value for this number. 

The zero temperature one-loop effective action is dehned then by 

IV = - In Z = -i In detL + 1cl(0) ^ 

= ~(y*a(ii))l..o (23) 

The term Cl( 0) hi /x^ = 04 ( 1 , L) In /x^ in the effective action W determines 
the one-loop beta function, this term describes renormalization of the one- 
loop logarithmic divergences appearing in the theory. 

2.3 Free energy for boson fields 

A hnite temperature held theory is dehned in Euclidean space, since for boson 
helds one has to impose periodic boundary conditions on imaginary time 
coordinate (antiperiodic boundary conditions for fermion helds respectively). 
A partition function is dehned by 

Z{p) = Tre-^^, (24) 

where if is a hamiltonian of the problem and (3 = h/T. Let us choose the 
lagrangian density p in the form 

d = + (25) 

where r is an imaginary time coordinate and L is a three dimensional spatial 
part of the density in the form (HD- The free energy of the system is dehned 
by 

FiP) = = -^\n{Np j D(j)exp{-dr j d^x (j)p(j))'^ , (26) 
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the integration is over all periodic fields satisfying ^(r + (3) = (j){T) (iVg is a 
normalization coefficient). As a result the eigenfunctions of p have the form 
exp(ira;,i)0A, where Un = 27rn//3 and L(l)\ = A0a- The free energy is thus 
equal to 




2(3 

' n=—oo A 

where we introduced ^-function 


/i" 


h d 




(27) 


+ 00 


C(^) — ^ + ■^) 


(28) 


n=—oo X 


and the parameter p with a mass dimensionality in order to make the argu¬ 
ment of the logarithm dimensionless (also see a previous section). 

Then it is convenient to use the formula 


C(^) = 


f + OO +°o 

j.s-1 


r(s) 


dtf 


EE- 

n=—oo A 




(29) 


and separate n = 0 and other terms in the sum. For n ^ 0 terms we 
substitute the heat kernel expansion for the operator L at small t 


e-"'=/<■(/; Li«) ~ at(/, L)(“ 


-3)/2 


(30) 


k=0 


and perform t integration, then we arrive at the high temperature expansion 
{(3 —>■ 0) for the free energy F\ 


F/n = -^cao) - Aa(0) i.i(,ia + (4^)=''= 


Os 


In 


ff3p\ 

04 / 

\27r) 

IOtt^ \ 


7 In 


27r 


E 

n>5 


(27r)3/2- 


/9' 


'4—n 


2 V 2 


n 


Cnin - 3) 


ao ai Ch(3) 02 1 

>490 /?3 47^3/2 p 2 24 


(31) 


Here = ak{I,L), Cr{s) = ^ Riemann zeta function, Ciis) = 

is a zeta function of an operator L, 7 is the Euler constant. The first 
two terms on the r.h.s. of m follow from the n = 0 term. 

The term 


TT^ ^, trJ vr^ 


= -V - 

(3^ 90 90 


( 32 ) 












is the leading high temperature contribution to the free energy. 

The classical limit terms due to the equality Cl( 0) = as can be rewritten 
as follows: 

T(-icU0)+a(0)ln5^)=T^l..^. (33) 

A 

The terms on the l.h.s. of (jSSI) yield a renormalized value of the terms on 
the r.h.s. of dSSl), since the sum on the righthandsight is generally divergent 
when the number of modes is inhnite. 

The term with 04 determines the part of the free energy that appears due 
to one-loop logarithmic divergences and thus it depends on the dimensional 
parameter /i as in the zero temperature case. 

2.4 Chiral anomaly in four dimensions 
for MIT bag boundary conditions 

Consider the Dirac operator on an n-dimensional Riemannian manifold 

5 = y (34) 

in external vector and axial vector helds. We suppose that and 
are anti-hermitian matrices in the space of some representation of the gauge 
group, is the spin-connection^. 

The Dirac operator transforms covariantly under inhnitesimal local gauge 
transformations (the local gauge transformation is D ^ exp(—A)D exp(A)): 

^\A^j, = [4^, A] 

hV, = d,\+[V^,\] 

D^D + [D,\] (35) 

and under inhnitesimal local chiral transformations (the local chiral trans¬ 
formation is D —> exp(z(p 75 )iA exp(i 9 ? 75 )): 

-f- [V)t, (f], 

D ^ D + z{D,j^ip}. (36) 

^The spin-connection must be included even on a flat manifold if the coordinates are 
not Cartesian. 
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The parameters A and tp are anti-hermitian matrices. 

First we adopt the zeta-function regularization and write the one-loop 
effective action for Dirac fermions at zero temperature as ^ 

W = -IndetD = -^IndetD^ = ^( 52 ( 0 ) + ^ ln(/i 2 )C 52 ( 0 ), (37) 

where 

CB.W = Tr(B-2-). (38) 

prime denotes differentiation with respect to s, and Tr is the functional trace. 
The following identity holds: 

a(s) = TrA-^ ^ 5a(s) = -sTr((524)24-^-i). (39) 

Due to the identity o 

5aC5^(s) = - (2sTr([D, A]D-2*-')) = -2s (Tt{[D-^% X])^ =0, (40) 

so the effective action (jSID is gauge invariant, 6\W = 0. 

The chiral anomaly is by definition equal to the variation of W under an 
inhnitesimal chiral transformation. Using o we obtain: 

ACdAs) = - (2!s'IV({B, 7y}B-"-‘)) = -to (iV)!"^®-"*)) . (41) 

and the anomaly reads 

^ = ^^^( 52 ( 0 ) = -2Tr(z7V5-'^)|.=o • (42) 

The heat kernel is related to the zeta function by the Mellin transforma¬ 
tion: 

poo 

Tr(i7V5"^') = r(s)-W . (43) 

Jo 

In particular, after the substitution of the heat kernel expansion m into 
the formula (HHl) we obtain 


A=-2an{t-f^^,D^). (44) 

The same expression for the anomaly follows also from the Fujikawa approach 

^ _ 

^The one-loop effective action is proportional to Planck constant h, in what following 
we put h= 1. 
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One can also derive the expression for the anomaly (I44|l from Schwinger’s 
effective action. One should start from an identity: 

/■+“ rif 

\nX = - -e-*^ (45) 

Jo t 

Then the change in the effective action due to chiral transformations can be 
written: 


•J 0 

/• + 00 _ /• + C30 f\ _ 

= -itr / { 7 V, 5 ^}e"*^' = 2itr7V / = 

Jo Jo 

= —2 liintri7®93 = —2an{i'y^(p,D‘^). (46) 

We impose local boundary conditions: 

= 0 , n_ = ^ (l - 7 ^ 7 „) , (47) 

which are nothing else than a Euclidean version of the MIT bag boundary 
conditions |2S|. For these boundary conditions II)! = n_, and the normal 
component of the fermion current '0^7nV’ vanishes on the boundary. Spectral 
properties of the Dirac operator for bag boundary conditions are intensively 
studied | 2 H 1 - 

Since D is a first order differential operator it was enough to fix the 
boundary conditions (gZl) on a half of the components. To proceed with a 
second order operator L = we need boundary conditions on the remaining 
components as well. They are defined by the consistency condition j29j : 

n-D'0|aM = O, (48) 


which is equivalent to the Robin boundary condition 

(V„ + R)n+V’|aM = 0, n+ = l(l + 7'57„) (49) 

with 

^ = -^n+L,,. (50) 

In the paper [H] the following expression for a coefficient 04 (Q, L) with 
an hermitian matrix valued function Q and conditions (0, Lab = 0 (flat 
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boundaries), S' = 0 was obtained: 


a^iQ, L) = ^(47r)-"/2| f {Q(60E./ + 60RE + 180E^ 

odU t Jm 

+ 5R^ - + 2R^,p^R>^'^P^)} 

+ f d"'~^xy/htT {Q{30E.n +^OxE-nX + ^OxE-n + QOE-riX 

JdM 

+ l8xX-.a^an + l2X-aRt‘anX + ^^^anXX’.a “ l2x^anX-.a 
+6[x^anX5 X:a\ + ^MX'.a, ^an] + 30[x, ^an:a] + + 30X-R;n} + 

+Q;n(-30E + 30x^X + 90x^ + 90Ex - 

-18X-.aX:a + 30X:aa " 6X:aX:aX + 30xi?) + 30xQ;/"}}. (51) 

To obtain the chiral anomaly in four dimensions^ with MIT bag boundary 
conditions one has to calculate the coefficient 04 (Q, L) (jBlj) with L = D^, 

Q = ix^(j) and substitute it into (m. We dehne - dy^ + [V^, W], 

= D^A^-D^A^, D^A^ = d^A„-TP^^Ap+[Vf,, A^]. The anomaly contains 
two contributions: 

A = Av + Ab. (52) 

In the volume part 

+60 [DpA,, 1/n - 60 DpD>^Dy^ + 120 {{DpA,, A"}, A^} 

+60 {DpA^^, Ay^} +120 Ayyy^ + 30 [[a^, a,], Ay 

+Wa {-45 i V^yp^ + 15 i A>^yp^ - 30 i {V^'^ApA^ + A^Ayp^) 
-120 i A^V^PA^ + 60 i A^A^pA^ - 60 [DyAjR''" + 30 {Dy>^)R 

-^ep,p^RP%0Rp^P^^ (53) 

only the DA — R terms seem to be new jH] (for flat space it can be found e.g. 
in ID). 

The boundary part 

= lon^o ^2 / d^x\/htT (12 i {Ab, ip}DaA^ 

180 (27r)^ Jqm V 

+24{ip, A^}{Aa, 414 - 60 4“, if] (14a - 4a, 4,]) 
+60{Dnif)DpAP^ (54) 

^In two dimensions (n = 2) the boundary part of the chiral anomaly with MIT bag 
boundary conditions is equal to zero [HI- 
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is new p. It has been derived under the two restrictions: S' = 0 and Lab = 0. 
Note, that in the present context, the hrst condition (S' = 0) actually follows 
from the second one {Lab = 0) due to (j50|l . 


3 Casimir effect for rectangular cavities 


3.1 Casimir energy of two perfectly conducting paral¬ 
lel plates 

The Casimir energy is usually dehned as 

E = (55) 

i 

where the sum is over all eigenfrequencies of the system. In what following 
we put h = 1. We start from the well known case of two perfectly con¬ 
ducting plates separated by a distance a from each other. In this case the 
eigenfrequencies oji are defined as follows: 


ojte = \J (vrn/a)^ + k1 + n = 1.. -1- oo (56) 

ujtm = \J{7m/ay + k^, n= 1.. + oo (57) 


^mainwave — 



(58) 


SO that the Casimir energy can be written as 



S is the surface of each plate. The hrst sum is equivalent to the sum over 
eigenfrequencies of the scalar held satisfying Dirichlet boundary conditions, 
the second sum is equivalent to the sum over eigenfrequencies of the scalar 
held satisfying Neumann boundary conditions. 

The expression for the Casimir energy written in this form is divergent. 
One has to regularize it somehow to obtain a hnite answer for the energy. 
Diherent methods were used for this purpose. In the present paper we suggest 
a method which makes calculations of determinants straightforward and easy 
to perform. 

By making use of an identity 





k 


(60) 


13 








we can see that up to an irrelevant constant the Casimir energy can be written 
in the form (we introduce a dimensional parameter /i by the same reasoning 
as in (j 2 ^ or pTjl h 



Now the expression for the Casimir energy is written in the standard Trln = 
InDet form, which is usual for one-loop effective actions in quantum held 
theory. The coefficient 04 ( 1 , L) is equal to zero for our current choice of the 
operator L and boundary geometry. 

At this point we introduce a regularization - we restrict integrations over 
momenta by some cut off K in the momentum space. The sums over n are 
also restricted as follows: 


+00 +00 -\-N +N 


n=l n=0 n=l n=0 


(62) 


The regularized Casimir energy is dehned by: 



It is convenient to perform a summation over n hrst. The following iden¬ 
tity holds: 


N N N 

^ln(?7,^ -I- {ka/irY) = ^ln(l -|- {ka^/ir'^n'^) + ^ln(n^). (64) 

n=l n=l n=l 

The hrst sum in (ini can be calculated in the N —> -|-cxo limit by use of an 
identity: 

+ (65) 

n=l 

The second sum in can be derived by use of a Stirling formula (which is 
exact in the large N limit): 

m ~ ( 66 ) 
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so in the large N limit it is possible to write: 

N 

= 2 In N\ = ln{2n) + f{N). (67) 

n=l 

In the large N limit the Dirichlet sum (El can be rewritten as: 

N 

ln(ri^ + {ka/irY) = ka + ln(l — exp(—2A;a)) — ln(A;a/7r) + f{N). (68) 

71 = 1 

The sum over Neumann modes can be rewritten as follows: 

N 

ln(?7,^ + {ka/irY) = ka + ln(l — exp(—2/ca)) + ln(/ca/7r) + f{N). (69) 

n=0 

It is possible to add any hnite number that does not depend on a to the 
regularized Casimir energy E^eg (1^1 (the force between the plates is being 
measured in experiments, so the energy can be dehned up to a constant). 
We add the surface term 

-7^ f dk e 2f(N) (70) 

(27r) Jo 

to the regularized Casimir energy Ereg (El- Doing so we obtain 

Sa S / \ 

Ereg = 2——^ / dk k^ + / d/c 2 lu(1 - exp(-2/ca)) (71) 

(27r) Jo [2'^) Jo ^ ' 

The hrst term in dHD is twice the regularized Casimir energy of the free 
scalar held since it can be rewritten as 

I- 

This term should be subtracted because we are interested in the change of 
the ground state energy when the plates are inserted into the free space. 

Next we perform the limits K +oo, N —> +oo. The Casimir energy is 

thus 

c r+oo o 2 

^ = ( 2 ^i = - 7 ^. (73) 

which is the well known result by Casimir ng. 

After elaborations we summarize the key points of the method, which is 
valid for the calculations in cylindrical cavities with arbitrary cross sections. 
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Suppose we want to calculate Trln of the second order operator + 

L^‘^\ where the dimensionalities of the operators are denoted by numbers. 
At zero temperature in our case of interest the operator describes a 
scalar field inside an infinite waveguide of an arbitrary cross section with 
Dirichlet or Neumann boundary conditions imposed. The eigenmodes of 
the operators and are denoted by and respectively. The 
following expression is finite (as can be seen from the heat kernel expansion): 


1 

2 





+ Af' 


,(i) 




To obtain the initial determinant one should add to (HH) the term 


(74) 


\ n ^ '®(^) + const (75) 

i 

The sum of (|7i|) and m generally has the following structure (to obtain 
the total Casimir energy the sum over indices k and p has to be performed): 


+ h‘Z’{\flP) + hZL + 2o(iV), (76) 

where 

— a convergent sum, (77) 

k 

the term dZZl) yields the energy of interaction for two flat parallel plates sep¬ 
arated by a distance a inside an infinite waveguide of the same cross section 
as these parallel plates (the walls of a perfectly conducting waveguide are 
perpendicular to two flat parallel perfectly conducting plates inside it). The 
term (HZl) yields an experimentally measurable contribution to the Casimir 
energy of the cavity (see Sec. 3.7 for details). 

The term /2 is equal to the self-energy of an inhnite waveguide 

when a —>■ cx). For rectangular cavities the term a\J~^/2 can be trans¬ 

formed to TrlnL^^^ in the same manner as in the beginning of this section 
(see a transition from (jnni) to (jnn) ). For the operator L 2 '’ 'w® repeat the step 
(EH) and continue this cycle until the first term in the righthandsight of (EH 
gets the form of the vacuum energy in an infinite space, i.e. the form EH- 
The terms describe the self-energies of two parallel plates inside 

the waveguide due to Dirichlet and Neumann modes, these self-energies do 
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not depend on a. For flat boundaries Dirichlet and Neumann boundary 
contributions to the Casimir energy cancel each other identically as can be 
seen from the Seeley coefficient ai Cl, for two parallel plates it can be seen 
from the expressions (inHD and (inni) . 

A contribution from the last term s{N) is proportional to (at 

zero temperature it is just the effective number of modes inside an infinite 
perfectly conducting waveguide, and thus it is not relevant to the energy of 
interaction between the two plates inside the waveguide) and 04 ( 1 , (this 
term is also not relevant when the interaction of the two plates inside an 
infinite waveguide is studied). 

To implement dZlD we used the following equality: 


^ (7rn/a)^ + 

J-J- (Tin/aY 

n=l ^ ' ' 


H-C50 

n(i+ 


n=l 




sinhay 


(78) 


3.2 Casimir energy of a perfectly conducting rectan¬ 
gular waveguide 

For a perfectly conducting rectangular waveguide the technical issues can be 
done in analogy with two parallel plates. We tacitly assume that the reader 
understood how the regularization is introduced in our method, so we will 
write only main steps without bothering too much on divergent form of some 
expressions. The Casimir energy for unit length is: 

+(™ 2 /fe)^+P?+P 2 ) (79) 

For TM modes ni and n 2 take positive integer values from 1 to +cxo, for TE 
modes ni and n 2 take positive integer values and one of them can be equal 
to zero (ui = ^2 = 0 corresponds to the main wave case). 

So the energy can be rewritten as: 

^ = i S /'/'%^ln(^(7mi/a)2 + (vms/fe)^+Pi+P 2 ) = 

= 4^ E //9^l>'((™.t‘ + (>^«^)“)+P?+P2). ( 80 ) 

m,712 = —00 ^ ' 

where t = h/a 01 t = a/b. Using formula for multiplication over ni we 
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obtain for the energy: 


H-C50 


E = - V 

Aab ^ 


n2=—oo 

H-cxd 


(27r) 


f f dp,dp2 




t 


Aab JJ (27r)^ 

712 = —OC 

+ 21n^l - exp(^-2^{T^n^jtY + 


2W (7rn2/t)2 + ^i^-(21n2) + 




P? +Pi 


( 81 ) 

(82) 

(83) 


A contribution from the term (—2 In 2) in (j82j) should be subtracted following 
the analysis of Section 3.1. The part (I83j) with the logarithm is hnite, it 
contributes to the hnite hnal answer for the Casimir energy. 

For the hrst term in (IH2D we get: 


f f dpidp2 

2ab ^ ]J (27r)2 

712 = —OO 


1 

2Eab 



■ = 


720i2o6 


(84) 

(85) 

( 86 ) 


because up to a numerical coefficient the expression (ISfill is just the same as 
the formula (pITI) . 

So the Casimir energy for unit length of a rectangular waveguide can be 
written as the sum of and (IHni) : 


d^waveguide(^^i b') 


720M 
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3.3 Casimir energy of a perfectly conducting rectan¬ 
gular cavity 

The Casimir energy in this case can be written as: 


E 




( 88 ) 


Using formula dZHD and technique described in previous subsections we ob¬ 
tain: 



The remaining terms should be calculated (using formula again) as 
follows: 


-I-C50 

E 

ni=—oo 




3/1 


TT 

4 ^ 


(90) 
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As a result for the Casimir energy of the cavity we obtain: 


Yj Yj Y1 

^cavity c) 43 ^ 486 48c (65 c)H~ 


-\-oo 


I E 


n2,n3=—00 


/ TTZ t 1 ~ 

'—00 


27r 


— 2 a 


7rn2 


+ 


- +P^ 

c / 


( 91 ) 


3.4 Relation to the argument principle 

An argument principle is a convenient method of summation over the eigen- 
modes of the system (see ra and |2I] for its applications). The argument 
principle states: 

i *^^^ 0 ) - (92) 

where ujq are zeroes and Uoo are poles of the function /(a;) inside the contour 
of integration. For the Casimir energy 0(a;) = uj/2. We choose 


/(^) 



■ 

2 sin 

-kl-kl 


-k'i-kl 


(93) 


in case of a scalar held satisfying Dirichlet boundary conditions on the plates. 
The contour lies on an imaginary axis, a contribution from the right semicircle 
with a large radius is negligible. A denominator is chosen in this form to 
remove uj‘^ = + ky from the roots of the equation /(a;) = 0. In this case 
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we proceed as follows: 



// 



A 

dvrz 


II 



2 sin — ky 


IJ 

S 


III 


I 



dk In — exp(—2fca) j + 


(94) 


(27r)2 


V 


I 


K 




Stt^ 


= -r + volume free space contribution+ 

1440a3 

+ surface contribution. 

Here ut = iCo. We see that the argument principle is in agreement with dZHI). 

3.5 Zeta function regularization for the cavity 

(■-function has already been discussed in this paper, so it is natural to describe 
regularization of the Casimir energy for the cavity in terms of (-function. 
Usually the Casimir energy is regularized as follows: 



(95) 


where s is large enough to make (jHHl) convergent. Then we should continue 
analytically © to the value s = — 1 , this procedure yields the renormalized 
hnite Casimir energy. In our case eigenfrequencies uji should be taken from 
(IHHI). So the regularized Casimir energy of the cavity Ecavityi.0,, b, c, s) can be 
written in terms of Epstein ^ 3 ( 4 ^ 4, U s) and Riemann (ij(s) zeta functions: 



( 96 ) 
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( 97 ) 


111 +00 / 


ni 


n2 


— + ^ + — 


ns 


-s/2 


ni,n2,n3=—oo 
+ 00 ^ 


Cnis) = 


n=l 




(98) 


The prime means that the term with all n, = 0 should be excluded from the 
sum. The reflection formulas for an analytical continuation of zeta functions 
exist: 


r(0>r-'yR(s) = r(-^)R<-‘>/y„(i - s) 

r( 0 ir“‘''^Zs(a, b, c; s) = (a 6 c)“‘r(i^^)ir**“‘')''^Z 3 3 - s) 


By use of reflection formulas o, (unnD one gets: 

abc 


../111 
Zsi-, Ti —11 — 
\a 0 c 


Cr(-i) = 


27r3 

1 

12 ' 


Z,{a,b,c-,4) 


The renormalized Casimir energy can therefore be written as: 

n / 7 \ nbc „ , , ,, TT /I 1 1\ 

b, c) - -Y^Z,(a, 6. c; 4) + -(- + - + -j, 


(99) 

( 100 ) 

( 101 ) 

( 102 ) 

(103) 


One can check that the formulas (UnnD and p[|) coincide identically and yield 
the Casimir energy for a perfectly conducting cavity. 


3.6 Geometric interpretation 

In this section we suggest a geometric interpretation of the main formulas in 
terms of geometric optics. This interpretation clarihes the physical meaning 
of the results (IHZI) , (ED obtained, which is always important for further 
generalizations in more complicated cases. 

Several geometric approaches - a semiclassical method [22], a worldline 
approach jSS] and a method of geometric optics jHlj have been introduced 
recently for the evaluation of the Casimir energies. Our formulas dHZD, (ED 
yield a simple geometric interpretation for the Casimir energy of the rectan¬ 
gular cavities in terms of geometric optics. 
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Optical contributions to the Green’s function of the scalar held with 
Dirichlet and Neumann boundary conditions have the form: 

1)’^a/A„(x, x') exp{iuJilni^, x')) (104) 

n 

Gopticaii^^^'^^i) = ^ V^n(x, x') exp(zo;iG(x, x')). (105) 


n 


Here G(x, x') is the length of the optical path that starts from x and arrives at 
x' after n rehections from the boundary. A„(x, x') is the enlargement factor 
of classical ray optics. For planar boundaries it is given by A„(x, x') = 1//^. 
From (HI it follows that for two parallel plates the Casimir energy of the 
electromagnetic held can be expanded as: 


E = 


S 


^H-oo 


dk In^l — exp(—2fca)^ = 


^H-OO 


dvr/c^ exp(—2an/c) 

(27r)^ 2an 


= -2aS 


dk 


+ CXD 


EE 

LOi n=l 


exp{—2anuji) 


2an 


(106) 


Here Yhwi ^ photon states (with frequencies uji = ^k^ + ky + k^) 

in an inhnite space. The righthandsight of (fTOOll can be written in terms of 
optical Green’s functions: 


H-OO 


EE 

n=l 


exp{—2anuji) 
2 an 


271 '^{G^pticaii^, X, a, iUi) + x, a, iui)) 

UJi 

(107) 


Note that terms with odd rehections from Dirichlet and Neumann Green’s 
functions cancel each other due to the factor (—1)*^ present in optical Dirichlet 
Green’s function. This is why only periodic paths with even number of 
rehections from the boundary l 2 n = 2,an enter into the expression for the 
Gasimir energy. 

Now consider the formula for the cavity (jOTj) (to obtain the Gasimir energy 
of a waveguide dHZl) in terms of optical Green’s functions the arguments 
are the same, just start from two parallel plates). Imagine that there is a 
waveguide with side lengths b and c. In order to obtain the rectangular cavity 
we have to insert two perfectly conducting plates with side lengths b and c 
(and a distance a apart) inside the waveguide. The eigenfrequencies that 
existed in a waveguide were equal to oowave = \/(^ 2 /^)^ + (ttus/c)^ + 
Only the photons with frequencies oJwave existed in a waveguide, and these 
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photons start interacting with the plates inserted inside a wavegnide. The 
optical contribntion to the Casimir energy arising from the interaction of 
these uJwave photons with inserted plates is equal to 


-271 {G^pUcali^, X, a, iuj^ave) + X, tt, iuj^^ave)) = 


71 


48a 


+ 



— exp 



(108) 


(this equality is obtained in analogy to (llObj) and ()l()7|l i. After comparison 
with dnn) it is straightforward to rewrite the Casimir energy of the rectangular 
cavity in terms of optical contributions : 


T'caiiitj/c) 

— 271 ^ ^ G^p^^^^i{'x,:x.,b,iujpiates))— (109) 

^plates 

- 271 E wave ) +Gopticaii^,^,a,i^^ wavs') ) • 

^wave 

Here we sum over all eigenfrequencies of the electromagnetic held in 3 cases: 
when there is an inhnite space (cui = \/kl + ky + two parallel plates 
{ujplates = y/(7m2/6)^ + au inhnite waveguide 

{i^wave = \/(7rn2/fe)2 + (vms/c)^ + {k^Y). 

The hrst two terms in (from may have the following geometric interpre¬ 
tation: from (innD it follows that 

^ = 7r^(Gj,ii^„;(x,x,6,i^mi«)+ G'^ncai(x,x,6,i^m«i)) (HO) 

^mw 

^ = TT ^(G^w(x, X, C, iuJrnw) + G^t*caz(x, X, C, iuj mtu)) ; (111) 

^mw 

where ujmw is an eigenfrequency of a main wave in a waveguide. So it is pos¬ 
sible to express Casimir energies of perfectly conducting rectangular cavities 
in terms of optical Green’s functions only. 

It is interesting that the Casimir energy of a perfectly conducting cavity 
can be written in terms of eigenfrequencies of the electromagnetic held in 
a free space, between two perfectly conducting plates and inside a perfectly 
conducting waveguide. 
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3.7 The experiment 

For the experimental check of the Casimir energy for the rectangular cav¬ 
ity one should measure the force somehow. We think about the following 
possibility: one should insert two parallel perfectly conducting plates inside 
an inhnite perfectly conducting waveguide and measure the force acting on 
one of the plates as it is being moved through the waveguide. The distance 
between the inserted plates is a. 

To calculate the force on each plate the following gedanken experiment is 
useful. Imagine that 4 parallel plates are inserted inside an inhnite waveguide 
and then 2 exterior plates are moved to spatial inhnities. This situation is 
exactly equivalent to 3 perfectly conducting cavities touching each other. 
From the energy of this system one has to subtract the Casimir energy of an 
inhnite waveguide, only then do we obtain the energy of interaction between 
the interior parallel plates, the one that can be measured in the proposed 
experiment (the subtraction of the term ()72|1 is just the same subtraction for 
two parallel plates). Doing so we obtain the attractive force on each interior 
plate inside the waveguide: 




dEattri,^) 6, c) 
da 


( 112 ) 


where 





TT 

Eattr 

(a, 6, c) = 

48a 

1 

+ O0 

r+oo 

1 

+ 4 

E 



n2,n3=—oo 



dp 


In 


1 — exp 


—2a 



(113) 


coincides with (Uni. We note that our formula (HH for the special case 
b = c coincides with the formula (6) in reference |HS1, there it was obtained 
using a diherent method and presented in a diherent mathematical form. 

To obtain the energy of interaction between the opposite sides of a single 
cavity one should subtract from the expression (ED) the Casimir energy of 
the same box without these two sides, i.e. one has to subtract from (inT|l the 
expression for the Casimir energy of a waveguide of a hnite length. To our 
knowledge the expression for the Casimir energy of a hnite length waveguide 
is not known up to now. 

It was often argued that the constant repulsive force (for a hxed cross 
section) derived from ET|l can be measured in experiment. However, without 
the subtraction just mentioned it is not possible to measure the forces in any 
realistic experiment, this is why it is not possible to use the expression (EH) 
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directly to calculate the force in the experiments. However, it can be used 
to derive a measurable in experiments expression for the force between the 
parallel plates inserted inside an infinite waveguide of the same cross section 
as the plates. 

Using the same technique as before it is possible to generalize our for¬ 
mulas (uni), (nn for the case of an infinite waveguide with an arbitrary 
cross section. The force between the two plates inside this waveguide can be 
immediately written : 


F{a) = 


da 
1 


Earbi^E) ^ ^ „ ln(l 6Xp( Q.a UJujave)') •) 


(114) 

(115) 


the sum here is over all TE and TM eigenfrequencies ujwave for the waveguide 
with an arbitrary cross section and an inhnite length. Thus it can be said 
that the exchange of photons with the eigenfrequencies of a waveguide between 
the inserted plates always yields the attractive force between the plates. 

To get the free energy Farb{ci, P) for bosons at nonzero temperatures /3 = 
1/T one has to make the substitutions (see Sec. 2.3, the formula (pTjl h 


p 


Pr 


2Tim 

~T 



1 


+ O0 

E 


m=—oo 


(116) 

(117) 


Thus the free energy describing the interaction of the two parallel perfectly 
conducting plates inside an infinite perfectly conducting waveguide of an 
arbitrary cross section has the form: 

+00 , _ 

Farbia, p) =-J2 

AfcD m=-oc 

+00 ^ _ 

E E + (118) 

^iNeum ^ 

where and Xf^eum eigenvalues of the two-dimensional Dirichlet and 
Neumann problems (a boundary here coincides with the boundary of each 
plate inside the waveguide): 

= -Xlofkix^y) 

fk{x,y)\aM = 0, 
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(119) 

( 120 ) 






( 121 ) 

( 122 ) 


A*^'9.(a:, V) = ->HNeum9i(x, y) 

Sgi(x,y) 


dn 


dM 


= 0 . 


The attractive force between the plates inside an infinite wavegnide of 
the same cross section at nonzero temperatnres is given by: 


F{a,P) = 


dFarb{a,(3) 


da 


1 




lotd 


--T 

R ^ 


utn 


(3 ^ exp(2au;rD) — 1 (3 ^ exp(2aci;Tv) — 1 


(123) 


ujtd 


u)tN 


Here ljtd = \/Pm + and ujtn = + A i 

The proof of these resnlts will be presented elsewhere. 


'iNeum' 
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